Few-boson dynamics in double wells: From single-atom to correlated-pair tunneling 



00 
O 

o 

(N 
C 



Sascha Zollner/'E] Hans-Dieter Meyer, 1 and Peter Schmelcher 1 ' 2 

1 Theoretische Chemie, Universitat Heidelberg, Im Neuenheimer Feld 229, 69120 Heidelberg, Germany 
2 Physikalisches Institut, Universitat Heidelberg, Philosophenweg 12, 69120 Heidelberg, Germany 

(Dated: November 22, 2007) 

We investigate few-boson tunneling in a one-dimensional double well, covering the full crossover 
from weak interactions to the fermionization limit of strong correlations. Based on exact quantum- 
dynamical calculations, it is found that the tunneling dynamics of two atoms evolves from Rabi 
oscillations to correlated pair tunneling as we increase the interaction strength. Near the fermioniza- 
tion limit, fragmented-pair tunneling is observed and analyzed in terms of the population imbalance 
and two-body correlations. For more atoms, the tunneling dynamics near fermionization is shown 
to be sensitive to both atom number and initial configuration. 
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The double well is a paradigm model for some of the 
most fundamental quantum effects, like interference or 
tunneling. Using ultracold atoms, it has become possible 
to study this system at an unprecedented level of control. 
This has lead, e.g., to the observation of Josephson os- 
cillations [l|, 0, 0] and nonlinear self-trapping [3, 0, 0] of 
Bose-Einstein condensates. In the first case, the weakly 
interacting atoms — prepared mostly in one well — simply 
tunnel back and forth between the two wells in analogy 
to a Josephson current. However, above a critical inter- 
action strength, the atoms essentially remain trapped in 
that well for the experimental lifetime even though they 
repel each other. On the few-body level, this resembles 
the situation of repulsive atom pairs, whose stability @ 
and dynamics have recently been observed. 

All of these effects are confined to the regime of rel- 
atively weak interactions, where the dynamics can be 
understood qualitatively (up to phases) by means of a 
single parameter: the number of atoms in one well. How- 
ever, interactions in ultracold atoms can be adjusted ex- 
perimentally over a wide range, e.g., via Feshbach reso- 
nances 0|. In particular, in one dimension (ID) one can 
tune the effective interaction strength at will by exploit- 
ing a confinement-induced resonance 0], which makes 
it possible to explore the limit of strong correlations. 
If the bosons repel each other infinitely strongly, they 
can be mapped to noninteracting fermions [10(] in the 
sense that the exclusion principle mimics the hard-core 
interaction. While local properties like the densities are 
shared with their fermionic counterparts, nonlocal as- 
pects such as their momentum distribution are very dif- 
ferent. Sparked also by its experimental demonstration 
111 . T3l. this fermionization has attracted broad interest 
(see [l3Lll3| and Refs. therein). 

In this Letter, we investigate the case where a few 
atoms are loaded into the same well and explore the 
tunneling dynamics as we vary the interaction strength 
from zero up to the fermionization limit. For two atoms, 
we show that the character of the tunneling changes 
from Rabi oscillations to correlated pair tunneling. Near 



fermionization, the strongly interacting atoms tunnel 
back and forth as a fragmented pair. For three or 
more atoms, the tunneling dynamics turns out to depend 
strongly on the atom number and the initial imbalance. 

Model and computational method The double-well 
dynamics is described by the many-body Hamiltonian 
H = Eti [M + U(xi)] + S E l<J S^Xi - x,). Here the 
double well U(x) = \x 2 + h8 w (x) is modeled as a su- 
perposition of a harmonic oscillator and a central bar- 
rier shaped as a Gaussian S w (x) = e~ x ' 2w /s/^irw (we 
choose w = 0.5 and h = 8, where harmonic-oscillator 
units are employed throughout.) The effective interac- 
tion resembles a ID contact potential but is mollified 
with a Gaussian <5 CT =o.05 so as to alleviate the well-known 
numerical difficulties of the S function. We focus on re- 
pulsive forces g £ [0, oo). 

Our goal is to investigate the few-atom quantum dy- 
namics in the crossover to the highly correlated fermion- 
ization limit g — ► oo in a numerically exact fashion. This 
is a challenging task, and most studies on the double- 
well dynamics so far have relied on two-mode models 
0, [H| valid for sufficiently weak coupling. Our ap- 
proach rests on the Multi-Configuration Time-Dependent 
Hartree method , a wave-packet dynamics tool which 
has been applied successfully to few-boson systems (see 
for details). 

From uncorrelated to pair tunneling To prepare the 
initial state \&(0) with a population imbalance — in our 
case, such that almost all atoms reside in the right-hand 
well — we make that side energetically favorable by adding 
a linear external potential — d ■ x (d > 0) and let the 
system relax to its ground state , I / g d>0 ' ) . For sufficiently 
large d, this amounts to preparing nearly all atoms in 
one well. To study their time evolution in the symmetric 
double well, in our simulations the asymmetry will be 
ramped down, d — » 0, within some time r > 0. 

Let us now study how the tunneling changes as we pass 
from uncorrelated tunneling (g = 0) to tunneling in the 
presence of correlations and finally to the fermionization 
limit (g — ► oo). It is natural to first look at the con- 
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Figure 1: (color online) Two-atom dynamics, (a) Relative 
population of the right-hand well over time, PB.(t), for differ- 
ent interaction strengths g — ( — ), g = 0.2 ( ), g = 4.7 

(■•■), and g = 25 (— • — ). (b) Snapshots of the one-body den- 
sity p(x) for different times t in the fermionized case g = 25. 
(All quantities in harmonic-oscillator units throughout, see 
text.) 

ceptually clearest situation where N = 2 atoms initially 
reside in the right-hand well. Absent any interactions, 
the atoms simply i?a6i-oscillate back and forth between 
both wells, which materializes in the percentage of atoms 
in the right well pr (t) = (0(a:)}^( t ) = J °° p(x;t)dx (p be- 
ing the one-body density) or, correspondingly, the pop- 
ulation imbalance 5 — pr — pl = 2pR — 1. By contrast, 
if the atoms repel each other, then the tunneling process 
will be modified, as can be seen in Fig.QJa). For g = 0.2, 
one sees that the tunneling oscillations have become a 
two-mode process: There is a fast (small-amplitude) os- 
cillation which modulates a much slower oscillation in 
which the atoms eventually tunnel completely (pr r* 0). 
In case g is increased further, we have found that the 
tunnel period becomes indeed so long that complete tun- 
neling is hard to observe. E.g., at g — 1.3 the period is as 
large as 2 x 10 3 . What remains is a very fast oscillation 
with only a minute amplitude - the two-body analogue of 
quantum self-trapping. As we go over to much stronger 
couplings (see g = 4.7), we find that the time evolution 
becomes more and more complex, even though this is 
barely captured in the reduced quantity pr [Fig. QIa)]. 
What is striking, though, is that near the fermionization 
limit (see g = 25) again a simple picture emerges: A fast, 
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Figure 2: (color online) Low-lying spectrum of two bosons 
in a double well as a function of the interaction strength g. 
Inset: Doublet formation with increasing g. 

larger-amplitude motion is superimposed on a slightly 
slower tunneling oscillation whose period roughly equals 
that of the Rabi oscillations. 

To get an understanding of the oscillations, Fig. [^ex- 
plores the evolution of the two-body spectrum {E m (g)} 
as g is varied. In the noninteracting case, the low-lying 
spectrum is given by distributing the N atoms over the 
lowest anti-/symmetric orbital of the trap. This yields 
the N + 1 energies {E m = Eq + mA' '}|J =0 , where 
A(°) = ei — cq is the energy gap between these two or- 
bitals, or the splitting of the lowest band. Assuming that 
for sufficiently small g still only N + 1 = 3 levels are pop- 
ulated, then the imbalance S(t) (and likewise pr) can be 
computed to be [ijj 

S(t) = S i cos(cjnii) + $12 cos(wi 2 i), (1) 

where u mn = E m - E n and 6 mn = 4(* m |9(a;)|* n )c m c„ 
are determined by the participating many-body eigen- 
states and their weight coefficients c m . At g = 0, due to 
the levels' equidistance, only a single mode with Rabi 
frequency cJoi = W12 = A^ ^ contributes. However, 
as the interaction is "switched on", the two upper lines 
El, 2 virtually glue to one another to form a doublet, 
whereas the gap to E increases (Fig. [21 inset) . For times 
t -C X12 = 2tt/uji2, we only see an oscillation with period 
Tqi <C T12, offset by 612, which on a longer timescale 
modulates the slower oscillation determined by U12. For 
small initial imbalances, | Co / C2 1 = | <^oi /<5i2 1 ^ 1; so for 
short times we observe the few-body analogue of Joseph- 
son tunneling. In our case of an almost complete im- 
balance, in turn, | <5i 2 1 dominates, which ultimately cor- 
responds to self-trapping, viz., extremely long tunneling 
times. These considerations convey a simple yet essen- 
tially exact picture for the two-body counterpart of the 
crossover from Rabi oscillations to self-trapping beyond 
the bare two- mode approach common for condensates 0. 

It is obvious that the two-frequency description above 
breaks down as the gap to higher-lying states melts, as 
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Figure 3: (color online) Top: Probability P2{t) of finding two 
atoms in the same well for g = 0,0.2,25. Bottom: Snap- 
shots of two-body correlation function P2(xi,x 2 ) at equilib- 
rium points, S(t) = 0, for g = (t = 44), g = 0.2 (t = 128), 
and g = 25 (t = 53) - from left to right. 



for g = 4.7. Concordantly, the dynamics becomes more 
complicated. However, in the fermionization limit (exem- 
plified for g = 25), the system becomes integrable again 
by mapping it to noninteracting fermions [lfj|. As an 
idealization, assume that at t — we put two (auxiliary) 
fermions in the ground state of the right well, where they 
would occupy the lowest two orbitals. Expressing this 
through the fermionic eigenstates |n)_ of the full sys- 
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where ll denotes occupation of the symmetric (a = 0) 
or antisymmetric (a = 1) orbital in band /?. Analyzing 
the corresponding energies, one finds that the frequen- 
cies contributing to the imbalance dynamics are exactly 
A' ' (the lowest-band Rabi frequency, corresponding to 
the longer tunneling period) and (the splitting of 

the upper band). This intriguing result states that only 
two modes determine the imbalance dynamics, so that 
the strongly repulsive atoms coherently tunnel back and 
forth almost like a single particle. As an illustration, 
snapshots of the density at different t are displayed in 
Fig. HUb). This demonstrates the tunneling of a frag- 
mented pair. 

In order to unveil the physical content behind the tun- 
neling dynamics, let us now investigate the two-body cor- 
relations. Noninteracting bosons simply tunnel indepen- 
dently, which is reflected in the two-body density (or cor- 
relation function) p 2 (xi,x 2 ). As a consequence, if both 
atoms start out in one well, then in the equilibrium point 
of the oscillation it will be as likely to find both atoms in 



the same well as in opposite ones. This is illustrated in 
Fig. [H which exposes p 2 at the equilibrium points and vi- 
sualizes the temporal evolution of the pair (or same-site) 
probability p 2 = L Xl . X2 > \ Pifeit x 2 )dx 1 dx 2 . As we in- 
troduce small correlations, the pair probability does not 
drop to 0.5 anymore - at g = 0.2 it notably oscillates 
about a value near 100%. This is apparent from the 
equilibrium-point snapshot of p 2 : Both atoms remain es- 
sentially in the same well in the course of tunneling. In 
other words, they tunnel as pairs. On top of this, Fig. [3] 
in hindsight also lays bare the nature of the fast (small- 
amplitude) modulations of Pn(t) encountered in Fig.QJa) 
by linking them to temporary reductions of the pair num- 
ber p 2 . Thus it is fair to interpret them as attempted 
one-body tunneling. As before, the time evolution be- 
comes more involved as the interaction energy is raised 
to the fermionization limit (cf. g — 25). The two-body 
correlation pattern is fully fragmented not only when the 
pair is captured in one well (corresponding, e.g., to the 
upper right corner Xi,x 2 > 0), but also when passing 
through the equilibrium point t = 53. Similarly, the evo- 
lution of p 2 {t) is governed by two modes, A^ ^ ± A^ 1 ), and 
over time p 2 passes through just about any value from 1 
(fragmented pair) to almost zero (complete isolation). 

Many-body effects Although having focused so far on 
the case of N = 2 atoms, the question of higher atom 
numbers is interesting from two perspectives. For one 
thing, it is fascinating because for g ^> 1 many results be- 
come explicitly A-dependent, including distinctions be- 



tween even/odd atom numbers [14]]. (The experimental 
preparation of definite N = 3, 4, . . . is feasible, if harder 
to achieve due to losses. In fact, the experimental setup 
in requires only an additional central barrier cre- 
ated by a Gaussian light sheet.) On the other hand, in 
a setup consisting of a whole array of ID traps like in 
[ill . fl2 |. number fluctuations may automatically admix 
states with N > 2. 



For N > 3, the weak- interaction behavior does not 
differ conceptually. In fact, Eq. J!]) carries over but 
with the sum now running over m < n < N . While 
the dynamics is no longer determined by strictly two fre- 
quencies, the separation of time scales (related to the 
formation of doublets in the spectrum) persists - ulti- 
mately, this should connect to the condensate dynam- 
ics valid for N 3> 1. Things become more intricate if 
we leave the two-mode regime, though. In particular, 
the fermionization limit reveals a clear N dependence 
(Fig. [4]) . Generally, an idealized state with N fermions 
initially in one well has contributions from all excitations 
1 1 £ o , . . . , llw-i )- ( a /3 = 0, 1 V/3) in the N lowest bands. 
Hence all tunnel split ting s A^) for each band are ex- 
pected to be present (l7| . Figure IDJa) conveys an im- 
pression of the complexity of the dynamics by exhibiting 
Pn(t) for N = 3,4. This somewhat erratic pattern may 
wash out the clear signature of the two-atom case upon 
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Figure 4: (color online) Many-body effects in the fermioniza- 
tion limit (g = 25). (a) Population of the right-hand well, 
PB,(t), for N = 3,4 atoms initially in one well. Bottom: Den- 
sity evolution p(x;t) for N—l = 2 (b) and N—l — 3 atoms (c) 
initially in the right-hand well if exactly one atom is present 
on the left. 



Finally, we mention that one may not only use the tilt d 
to load the atoms into one well, but also to study tunnel- 
ing oscillations in asymmetric wells in order to actively 
tune the tunneling. A detailed investigation [l3| reveals 
that, for medium g, single-particle tunneling can be res- 
onantly enhanced if the right well is lowered enough to 
compensate the interaction-energy shift. In the fermion- 
ization limit, in turn, single-atom tunneling turns out 
resonant already for d = 0, while tuning d makes other 
resonances accessible. 

In conclusion, we have performed an ah initio investi- 
gation of the full crossover from uncorrelated to fermion- 
ized tunneling of a boson pair in a double well. Remark- 
able features of this pathway are the strongly delayed 
pair tunneling encountered for medium interactions and, 
in the fermionization limit, fragmented pair tunneling at 
the Rabi frequency. Having pushed the notion of tunnel- 
ing toward strongly interacting systems, this opens up 
intriguing perspectives, ranging from resonantly tuning 
the tunneling to considering multi-well setups. 

Financial support from the Landesstiftung Baden- 
Wiirttemberg is acknowledged by P.S. and S.Z. The au- 
thors thank S. Jochim and 0. Alon for discussions. 



averaging over an array. In an experiment, it is therefore 
desirable to reduce number fluctuations, e.g., by having 
sufficiently high barriers in between different copies of 
the double well. 

In the context of many-body effects, it is interesting to 
consider what happens if not all N > 3 atoms are pre- 
pared in one well, but rather, say, N — 1 in one well and 
one in the other. Paraphrased in the case N = 3, this is 
the question of the fate of an atom pair if the target site 
is already occupied by an atom. The striking answer, as 
evidenced in Fig.[4][b), is that the process can be viewed 
as single-atom tunneling on the background of the sym- 
metric two-atom ground state. The tunneling frequency 
in the fermionization limit is simply the tunnel splitting 
A^ 1 ) m 27r/40. This has the intuitive interpretation of a 
fermion which — lifted to the band (3 = 1 — tunnels inde- 
pendently of the two lowest-band fermions. From that 
point of view, it should come as no surprise that adding 
another particle destroys that simple picture. In fact, 
Fig.SKc) reveals that if we start with N — 1 = 3 atoms on 
the right, then the tunneling oscillations appear erratic 
at first glance, and a configuration with three atoms per 
site becomes an elusive event. (E.g., at t sa 22, three 
atoms are on the left site, whereas at t « 44, 72 three 
atoms are on the right.) In the spirit of the Fermi map 
above, this can be understood as superimposed tunneling 
of one atom in the first excited band (A^) and another 
in the second band (A^ 2 ) w 27r/15), while the remaining 
zeroth-band fermions stay inactive. 
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